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Abstract. The variable flux of TeV gamma-rays detected 
from Mkn 421 and Mkn 501 requires the presence of high 
energy electrons, which could in principle produce large 
numbers of electron/positron pairs, leading to an electro- 
magnetic cascade. We point out that this scenario can be 
avoided if electrons are accelerated to high energy rectilin- 
early, rather than being injected isotropically into a blob, 
as in most of the models of the GeV gamma-ray emis- 
sion. By balancing linear acceleration by an electric field 
against inverse Compton losses in the radiation field of the 
accretion disk we calculate the emitted spectra and find 
the conditions which must be fulfilled in order to exclude 
the development of electromagnetic cascades during ac- 
celeration. Assuming these to be fulfilled, we show that 
the maximum possible photon energy is approximately 

2/5 

lOMg TeV, where M% is the mass of the central black 
hole in units of 10 s M Q . In addition we compute the op- 
tical depth to absorption of TeV photons on a possible 
isotropic scattered component and on the observed non- 
thermal radiation (in the case of Mkn 421) and find that 
TeV photons can escape provided the nonthermal X-rays 
originate in a jet moving with a Lorentz factor 7b 8. 
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1. Introduction 

Gamma-rays of energy up to 10 GeV have been ob- 
served from many active galactic nuclei by the EGRET 
experiment on the Compton Gamma-Ray Observatory 



(von Montigny et al. 1995). Although hadronic models 
have been proposed (Mannheim & Biermann 1992) most 
models assume the gamma-rays are produced by inverse 



Compton scattering off electrons contained in a blob of 
plasma which itself is in relativistic motion in the general 
direction of the observer (Maraschi et al 1992; Dermer 
& Schlickeiser |l993| ; Blandford & Levinson yj395|) . These 
models differ in their assumptions about the soft target 
photons (Schlickeiser 1996), but they all assume the elec- 
trons to be isotropically distributed in the rest frame of 
the blob. The gamma-ray luminosity they predict is rel- 
ativistically boosted, so that the observer receives a flux 
which exceeds that measured in the rest frame of the blob 
by several orders of magnitude. 

In the TeV regime two blazars - Mkn 421 and Mkn 501 
- have been observed to emit gamma-rays of energy 



>4TeV (Punch et al. |l992j , Quinn et al. |1995| ). Of those, 
Mkn 421 is also an EGRET source and its spectrum is con- 
sistent with an unbroken power law spectrum of photon 
index ~ —2 between GeV and TeV energies, which makes 
it unlikely that these photons have been reprocessed by 
an electromagnetic cascade. Electrons of Lorentz factor 
exceeding 10 6 are required in the blob frame in order to 
explain the observations, and, if isotropically distributed, 
these would inevitably undergo photon scattering events 
well into the Klein-Nishina regime. Such events produce 
pairs, either directly via the triplet pair production process 
(Mastichiadis 1991) or indirectly via photon-photon pair 
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production. Thus, although there are no detailed calcula- 
tions of the development of a cascade in the anisotropic 
radiation field seen by a relativistic blob, a problem is 
quite likely to arise. 

In view of this, we present here a model in which elec- 
trons are accelerated to high energies rectilinearly in an 
electric field. For certain values of the potential, acceler- 
ation can be balanced by inverse Compton scattering on 
the accretion disk photons which occurs in the Thomson 
regime, thus avoiding the problem of a cascade. Since we 
are interested in the outgoing gamma radiation, we follow 
electrons which propagate outwards and calculate their 
energy as well as that of the associated photons as a func- 
tion of height above the accretion disk. Adopting a simple 
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model for the disk radiation we find that close to the disk, 
Thomson losses limit the particle energy to a relatively 
small value, whereas far from it the losses are smaller and 
particles reach higher energies before their acceleration is 
saturated. Lower energy gamma rays thus originate close 
to the disk whereas photons of TeV energy come from a 
distance of about 1000 gravitational radii away from the 
black hole. 

The system of a jet emerging along the rotation axis 
of an accretion disk which underlies our calculations is de- 
scribed in Sect. ||. The saturation of electron acceleration 
by inverse Compton scattering and the resulting spectrum 
is computed in Sect. ||. In the same Section we give also 
the necessary conditions which make the inverse Comp- 
ton scattering to occur in the Thomson regime. Finally 
we discuss the implications of our model in Sect. [| 

2. The disk/jet geometry 

We adopt a simple one-dimensional model of a jet in which 
there are localised regions where particles can be acceler- 
ated by an electric field. There are several ways of real- 
ising such a configuration. Magnetic reconnection (e.g., 
Schindler et al 1991 , Haswell et al 1992| ) leads to potential 
drops aligned along the magnetic field, but is amenable 
only to rough estimates. Almost rectilinear acceleration 
can occur for some trajectories as a result of an encounter 
with a relativistic oblique shock front, such as may be ex- 
pected within a relativistic jet (Begelman & Kirk 1990| ) 
and is an integral part of the 'surfatron' process (Kat- 
souleas & Dawson 1983) which has as yet not been ap- 
plied in an astrophysical context. Alternatively, an elec- 
tric field component along the jet axis may be induced 
close to the black hole (Blandford & Znajek |1977|, Mac- 



donald & Thorne[l98|, Bednarek & Kirk |1995 ) or arise m 
the sheath between a relativistic jet and its surroundings 
(Bisnovaty-Kogan & Lovelace 1995| ). 

The orientation of accelerating regions is important 
for the computation of the losses suffered by the particles. 
We will assume it to be along the jet axis. The maxi- 
mum potential which can be realised within an accelera- 
tion region scales with the local magnetic field strength 
and the coherence length scale of the region. For a con- 
stant opening angle of the jet, the dominant magnetic 
field component falls off inversely as the distance zr g 
from the black hole. Here r g is the gravitational radius 
r g = 2GA/bh/c 2 of a black hole of mass A/bh- Li this 
case, the coherence length is likely to be proportional to 
z, so that the potential available in an acceleration zone is 



independent of position in the jet (Bednarek et al 1995). In 



this paper, we treat a more general situation in which the 
maximum potential AV varies inversely as some power 
a of z, i.e., AV = Voz~ a . The coherence length will 
be taken to be a fixed fraction (< 1) of the distance 
from the black hole, so that the electric field is given by 
E(z) = Vqz~°'~ 1 / (r g {3). Although it is to be expected that 



the jet moves relativistically with bulk Lorentz factor ~ 10 
we assume the acceleration zones to be stationary in the 
'lab.' frame, as would be expected if they are associated 
with shock fronts caused by an obstacle at the edge of the 
jet. The particle dynamics are then described simply in 
terms of the electric field, and have no superimposed bulk 
motion until they have emerged from an acceleration zone 
and have had enough time to isotropise. 

The radiation fields which may be present and lead 
to inverse Compton losses are of two different kinds - 
an ambient almost isotropic field arising from backscat- 
tering of the luminos ity of the entire black hole/di sk/jet 
system (Sikora et al 1994 , Blandford & Levinson [1995 ) 
and an anisotropic component of photons which come di- 
rectly from the disk. Objects which emit TeV photons 
must have relatively weak isotropic radiation fields, so 
that in this paper we will concentrate on the direct radia- 
tion from the disk (see Sect. |]). An important property of 
an accretion disk is that its surface temperature falls off 
with increasing radius r. To include this effect we adopt a 
Shakura/Sunyaev profile and assume the disk emits black- 
body radiation at temperature T(r) — Tj n (ri n /r) 3 / 4 where 
Ti n is the temperature at the inner edge of the disk as- 
sumed to lie at r = r ln = 3r g . Thus, two parameters are 
required to specify the disk, and these we choose to be 
the luminosity ^Edd of the disk in units of the Eddington 
luminosity: £ Edd = 1.5 x lO" 35 ^^ 
the black hole, expressed in units of 10 s M 



1.5 x 10- 35 T ; t7v and the mass M 8 of 



3. Inverse Compton scattering 

The Lorentz factor of electrons accelerated rectilinearly by 
the electric field described in Sect. 2, is limited under cer- 
tain conditions by Compton scattering. The radiation re- 
sponsible for this can in principle either come directly from 
the disk, or have been scattered by surrounding material 
into a quasi-isotropic distribution. The Thomson losses of 
a particle moving along the jet axis in the anisotropic ra- 
diation field of the disk are straightforward to compute 
e.g., by generalising the results of Protheroe et al ( |l992| ): 



7= 0.25 x 7 2 ^ Edd M s - V 



3 s^ 1 



(1) 



Using this formula, one finds that losses balance gains in 
the electric field at an equilibrium Lorentz factor given by 



7cq (z) = 8.82 x lO-VV^r 172 ^ 2 



(2) 



where Vq is measured in volts. Losses on the isotropic radi- 
ation field depend on the radius typical of the scattering 
material R sc and on its optical depth r sc . We find that 
disk radiation dominates the energy losses for distances 



z < 80oo(iWipc) 2/3 (Wio~ 3 r 1/3 ^ 8 ~ 2/3 



(3) 



A limit on the column density of matter in Mkn 421 can be 
obtained from X-ray observations: rtn < 1.45 x 10 20 cm -2 
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(Fink ct al 1991) corresponding to a Thomson optical 
depth r sc < 1CP 4 . According to Eq. (3) this means that 
losses are dominated by the radiation field of the disk. 

Equations (1) and (2) are valid provided the electrons 
scatter in the Thomson regime. Approximating the black- 
body spectrum by monoenergetic photons of energy Zk^T ', 
this condition implies the requirement 



3k B T in z~ 



3 / 4 {(3z/r) 3 / 4 [l-(l + (z/r) 2 )- 1 / 2 ]} 7 <mc 2 (4) ^ 



for photons from a point at radius r on the disk. The quan- 
tity in braces on the left-hand side reaches a maximum of 
roughly 0.7 at r 2z, leading to 



7 < 1.4 x 10 4 (M 8 /W/ 4 ^/ 4 



(5) 



In order to channel energy efficiently into gamma-rays, 
the potential must accelerate electrons to an energy high 
enough for losses to be important. If this is not the case, 
the acceleration process puts energy primarily into parti- 
cles, which may subsequently isotropise and radiate, as in 
other models. Thus, the condition 7 cq mc 2 = eAV sepa- 
rates the region in which energy is put directly into radi- 
ation, which we call the 'radiation dominated' zone, from 
the 'particle dominated zone' where energy is first trans- 
ferred by the acceleration mechanism into particles, which 
cool only after acceleration has ceased. 

Figure 1 shows the constraints on the parameter space 
for black- hole mass Mg = 1, disk luminosity ^Edd = 0.1 
and f3 = .1. In addition, the dependence of 7 cq on z is 
shown for various values of the potential Vq and for a — 0. 

If a particle is accelerated linearly at a point in the jet 
within the permitted region lying between the dashed and 
solid lines in Fig. 1, its equilibrium Lorentz factor is given 
by Eq. (2). If, however, the electric field is too strong, and 
the particle finds itself above the Klein- Nishina boundary, 
the photons it produces will start a pair avalanche sim- 
ilar to that described by Bednarek & Kirk ( 1995 ). As a 
result, we can expect that the particle will not achieve the 
Lorentz factor given by Eq. (2). On the other hand, a weak 
electric field places the particle in the particle dominated 
region. Here the total available potential in an acceleration 
region is insufficient to produce particles of Lorentz fac- 
tor given by Eq. (2). Thus, the maximum possible Lorentz 
factor 7 max is achieved for jet parameters within the per- 
mitted region and at the highest z compatible with this 
restriction. This value depends only weakly on the mass 

18^-2/5,1/5 
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Fig. 1. The parameter space for rectilinear acceleration. The 
solid line shows the constraint 7cq = eAV/ (mc 2 ) dividing radi- 
ation and particle dominated zones. The dashed line shows the 
boundary between the Thomson and Klein-Nishina scattering 
zones. These boundaries are independent of the dependence of 
the potential along the jet and are plotted for j3 = 0.1. The 
dotted lines show 7 cq for Vb = 10 11 and Vo = 5 X 10 12 V. These 
lines have the slope 1 — a/2 and are shown for AV = constant, 
corresponding to a = 0. 



approximately given by r = r g z. Thus, e s oc z~ 3 / 4 and 
e av oc z~" +5 ' 4 . The energy converted into gamma-rays by 
a single particle passing through a distance dz is equal to 
the energy extracted from the electric field. If dn photons 
are thereby produced, one can write e av dn = ecE(z)dz oc 



l dz. As shown by Protheroe et al (1992) the photon 



spectrum dn/de can be found by convolving dn/de av with 
the real distribution p(e,e av ) of emitted photons. For the 
cases which are of interest here, however, this does not 
affect the spectral slope, which is given by 



(dn/de) oc 



-(4a-10)/(4a-5) 



(6) 



If we assume equal numbers of particles are present in the 
acceleration zones over a range of heights z\ < z < Z2, 
then we obtain a power-law spectrum between the photon 
energies e\ = e(z\) and ti — t(zi) where 



of the black hole 7max = 1.1 x WM^t^^ and is <z) „ Lfl x m -7y Q ^ M -Vy i^-a+5/4 eV 



achieved at a height z max = 1.6 x 10 5 M^ /5 £^ d j3 4 / 5 . 

An analytic expression can be found for the radiated 
photon spectra in the "radiation dominated" zone using 
the method of Protheroe et al ( |1992|) . The average energy 
e av of a photon produced at height z is approximately 
7 2 q e s , where e s is the energy of the dominant soft pho- 
tons, which are those originating at a point on the disk 



(7) 



For a — 0, a spectrum of index —2 is obtained. The max- 
imum possible photon energy is then 



3.4 x 10 13 ^ 3/5 4dd^ 8 2/5 eV ■ 



(8) 
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4. Discussion 

In our model, the gamma radiation observed from blazars 
is produced by electrons subject to rectilinear accelera- 
tion which is balanced by losses due to inverse Comp- 
ton scattering on accretion disk photons. In this case, the 
TeV 7-ray emission from an object such as Mkn 421 orig- 
inates in the upper part of the radiation dominated zone 
of the jet at a distance of about 1000 gravitational radii 
from the black hole. If the disk radiation dominates over 
that of the isotropic component at even larger distances 
(z « 6000), the maximum achievable photon energy for 
P « ^Edd ~ 0.1 is about 10 TeV. 

There are two interesting features of the proposed 
model. Firstly, the inclusion of an explicit acceleration 
mechanism leads naturally to a prediction of the site at 
which high energy gamma rays are produced. Secondly, 
the fact that the electron-soft photon collisions occur in 
the Thomson regime guarantees that the resulting gamma 
rays are emitted below the threshold for pair production, 
thus avoiding the complications of an electromagnetic cas- 
cade. 

However, severe additional constraints must be ful- 
filled if these photons are to leave the source and avoid 
absorption on other ambient radiation fields. Two types 
of field are important: (i) an isotropic component arising 
from scattered disk radiation and (ii) the observed non- 
thermal radiation. Of these, (i) is not directly observed in 
the blazars of relevance here. Thus, there is little informa- 
tion concerning how much matter surrounds the accretion 
disk (such as column density or typical distance scale). 
Detailed computations of the optical depth for 7-ray pho- 
tons escaping from the accretion disk in Mkn 421 have 
been performed by Bottcher & Dermer (1995). However, 
these authors assume a Thomson optical depth around the 
disk of 0.02 out to a distance of 0.1 pc. This is consistent 
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with ROSAT observations (Fink et al |1991j ) only if the X- 
rays are produced at zr g > 0.1 pc. If, as we assume, this is 
not the case, r sc < 10~ 4 and the absorption of TeV 7-rays 
by the isotropic component (i) is negligible. 

In case (ii), the nonthermal radiation is directly ob- 
served. According to our model, these photons are pro- 
duced in the particle dominated zone at z > 1000 where 
the electrons have time to isotropise and form a relativis- 
tic radiating blob as in other models of the GeV gamma- 
ray emission. The optical depth to absorption of the TeV 
gamma-rays in such a blob is given by 



16J 2 

T(£) = ^ 



de b J d(cosfl b )F(2 7b e b ) 



x (1 - cos6> b )er 77 [e/(27 b ),e b , 



(9) 



where D is the distance to Mkn 421, i? b is the radius of the 
blob, F(e) is the observed photon flux at energy e [which 
we approximate as a power law with index —2 below and 
—2.4 above 10~ 4 MeV, with a normalisation taken from 



,1 1 1 



i — 1 — 1 — 1 — 1 — r 

Mrk 421 




40TeV 



15.5 16 16.5 17 17.5 

log R. [cm] 



Fig. 2. The requirements for escape of 4 TeV and 40 TeV 
7-rays through the nonthermal radiation of a relativistic blob 
in Mkn 421. The Lorentz factor 75 and radius ii;, of a radi- 
ating blob above which the optical depth for 4 TeV (full line) 
and 40 TeV (dotted line) 7-rays is less than unity. The dashed 
line shows the minimum 7b which permits an X-ray variability 
time scale of one day. 

er 77 (e, e',0) is the cross section for e pair production of 
a photon of energy e off one of energy e' with 6 the angle 
between the directions of travel (Jauch 



Rohrlich |1955D . 
The full line in Fig. 2 shows the parameters of the blob 
in Mkn 421 (7b, Rb) for which the optical depth for 4 TeV 
7-rays is equal to unity. An additional constraint on the 
Lorentz factor of the jet arises from the X-ra y variability 
time scale t val ss 1 day (e.g. Fink et al 199l| ). Assuming 
that the blob moves directly towards the observer we find 



7b > i? b /ctvar > 3.8 x 10~ 16 i? b [cm] 



(10) 



the observations (Fink et al 1991, Makino et al 1987)] and 



which is shown in Fig. 2 by the dashed line. The region 
in Fig. 2 above the full and dashed lines gives the per- 
mitted parameters of the blob in Mkn 421 for which the 
7-rays with maximal energy of 4 TeV can escape. From 
it we may also constrain the Lorentz factor of the blob: 
7b > 7.5. For reasonable values of 7b (of the order of tens) 
and an opening angle of the jet of the order of a few de- 
grees (which determines i?b), the blob is located inside the 
particle dominated zone of the jet, in agreement with our 
assumptions. 
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